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COMPARATIVE  STUDY  OF  THE  AXIAL  AND  AZIMUTHAL  BUNCHING 
MECHANISMS  IN  ELECTROMAGNETIC  CYCLOTRON  INSTABILITIES 

! 

I.  INTRODUCTION 

12  3 

In  the  late  1950s,  Twiss,  Schneider,  and  Gaponov  independently 

proposed  a radiation  mechanism  involving  electrons  gyrating  in  an  external 

magnetic  field.  This  phenomenon  known  as  the  electron  cyclotron  maser 

4-29 

instability  has  been  the  subject  of  continuing  activity-  It  is  the 

relativistic  mass  dependence  of  the  electron  cyclotron  frequency  which 

causes  electrons  to  bunch  azimuthal ly  in  their  gyration  orbits,  and 

consequently  to  drive  the  instability.  Recently,  intense  research 

activities  to  employ  this  interaction  has  occurred  in  the  development  of 

a new  microwave  source  called  the  gyrotron  or  electron  cyclotron  maser. 

23 

Results  reported  so  far  are  impressive  in  terms  of  peak  power 
(1  GW  at  X = 4 cm)  and  in  terms  of  efficient  production  of  large  c.w. 
power*®’*®  (12  kW  at  X = 2.8  mm  and  1.5  kW  at  X = 0.9  mm).  The  high 
power,  high  efficiency,  short  wavelength  radiation  produced  by  gyrotrons 
can  be  employed,  among  other  applications,  to  achieve  local  or  global 
heating  of  tokamak  plasmas.3®’3* 

• 32  33-41 

On  the  other  hand,  Weibel  and  later  authors  have 

investigated  a different  electromagnetic  instability  driven  by  the 

anisotropic  velocity  distribution  of  electrons.  In  contrast  to  the 

cyclotron  maser  instability,  this  instability  is  nonrelativistic  in 

Note:  Manuscript  submitted  September  1,  1977. 
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nature  and  will  be  referred  to  as  the  Wei bel -type  instability  in  the 

present  paper.  In  the  absence  of  wave-particle  resonances,  the  basic 

42 

mechanism  of  Weibel-type  instability  is  axial  electron  bunching  in  the 
direction  of  wave  propagation,  caused  by  the  ^ x Lorentz  force, 
where  _v  is  the  electron  velocity  component  perpendicular  to  the  wave 
vector  and  is  the  wave  magnetic  field.  In  the  presence  of  an  electron 
thermal  spread,  the  distribution  of  axial  velocity  will  tend  to  spoil 
the  phase  synchrony  discussed,  thereby  reducing  the  growth. 

There  are  fundamental  differences  between  the  two  instabilities. 

In  addition  to  the  qualitative  difference  in  bunching  mechanisms,  the 
Weibel-type  instability  does  not  even  require  an  external  magnetic  field, 
except  to  define  an  axis  for  the  velocity  space  anisotropy.  Yet,  there 
are  also  similarities  between  the  two  instabilities.  Both  will  take 
place  for  right-hand  circularly  polarized  wave  propagating  along  the 
external  magnetic  field  and  the  unstable  spectra  both  scale  with  the 
electron  cyclotron  frequency.  From  such  similarities,  one  expects  a 
close  relationship  between  the  two  instabilities.  In  fact,  there  have 
been  speculations  among  gyrotron  researchers  that  the  Weibel-type 
instability  might  also  be  exploited  for  microwave  generation,  or  may 

4 

even  be  responsible  for  some  observations  in  gyrotron  experiments. 

From  the  literature  on  the  two  types  of  instabilities,  one  finds  that 
investigations  of  these  two  subjects  have  evolved  separately.  Most 
studies  concerning  the  Weibel-type  instability  have  been  based  on  non- 
relativistic  models  and  thus  do  not  include  the  cyclotron  maser  instability. 
On  the  other  hand,  most  treatments  on  the  cyclotron  maser  instability  have 
concentrated  on  the  azimuthal  bunching  mechanism.  Hirshfield,  et  al6 
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briefly  discussed  the  relative  importance  of  the  terms  originating  from 
the  two  mechanisms  in  their  result.  They  found  that  the  azimuthal 
bunching  mechanism  dominates  when  the  wave  phase  velocity  exceeds  the 
speed  of  light  and  vice  versa,  a conclusion  we  will  verify  from  a 
different  point  of  view  and  elaborate  further  in  the  present  treatment. 

This  study  is  motivated  by  recent  interest  in  the  cyclotron 
maser  instability.  Its  purpose  is  to  examine  the  two  bunching  mechanisms 
under  a unified  physical  picture  and  thereby  clarify  the  physical 
relation  between  the  two  types  of  instabilities.  Our  analysis  will  be 
based  on  the  following  electron  distribution  function  commonly  adopted 
for  cyclotron  maser  studies, 

fo  = 6(pi'p±o)  6(pz)/2itpl  (1) 

where  Pj_»Pz  are  the  transverse  and  axial  momenta,  respectively,  p^Q  is 
a constant  and  S(x)  is  the  Dirac  delta  function.  Eq.  (1)  is  a realistic 
beam  frame  representation  of  the  electron  beams  used  in  most  gyrotron 
experiments.  For  such  a cold  distribution  function,  wave-particle 
resonances  are  absent.  Our  results  are  thus  strictly  valid  for 
sufficiently  monoenergetic  electron  distribution  functions  which  can 
be  approximated  by  Eq.  (1). 

We  find  from  the  simple  model  to  be  presented  that  there  is  indeed 
a strong  coupling  between  the  bunching  mechanisms  associated  with  the  two 
instabilities.  They  are  so  related  and  inseparable  that  only  a unified 
physical  interpretation  can  fully  describe  the  physical  processes  involved 
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in  either  instability.  Except  for  the  special  cases  of  zero  external 
magnetic  field  or  infinite  wavelength,  the  two  mechanisms  are  always 
simultaneously  present  and  compete  rather  than  reinforce  one  another. 
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By  comparing  the  two  mechanisms,  we  obtain  a criterion  for 
determining  the  relative  importance  of  the  two  mechanisms  and  identify 
the  two  types  of  instabilities  as  the  fast-wave  and  slow-wave  branches 
of  the  right-hand  circularly  polarized  wave  without  actually  solving 
the  dispersion  relation.  From  the  comparison,  it  can  also  be  seen  that 
the  two  mechanisms  interact  in  such  a way  that,  for  a given  k2>  only  one 
branch  can  be  unstable.  Another  significant  feature  resulting  from  the 
comparison  is  that  the  (relativistic)  azimuthal  bunching  mechanism  may 
domi nate  at  "nonrealativistic"  electron  energies,  while  the  (nonrelativistic) 
axial  bunching  mechanism  may  dominate  at  "relativistic"  electron  energies. 

In  Section  II  we  compare  in  detail  the  azimuthal  and  axial  bunching 
mechanisms.  In  Section  III  we  illustrate  the  conclusions  of  Section  II  with 
numerical  examples.  Section  IV  gives  further  discussions. 
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II.  COMPARISON  BETWEEN  AXIAL  AND  AZIMUTHAL  BUNCHING  MECHANISMS 


Consider  an  infinite  and  uniform  plasma  (density  n)  immersed 
in  a uniform  magnetic  field  (B^).  For  waves  propagating  along  the 
magnetic  field,  it  is  well  known  that  the  electromagnetic  and  electro- 
static effects  are  decoupled,  resulting  in  purely  electromagnetic 
(7-E_  = 0)  or  electrostatic  (VxJE  = 0)  modes.  The  mode  of  interest  here 
is  the  electromagnetic  mode  with  right-hand  polarization  (i.e.,  electric 
field  rotates  in  the  same  sense  as  electron  gyration)  and  with  wave 
requency  (or  growth  rate)  comparable  to  the  electron  cyclotron  fre- 
quency (or  plasma  frequency).  In  such  a case,  the  plasma  ions  simply 
provide  a charge  neutralizing  background  and  their  dynamics  can  be 
neglected.  For  the  case  of  a pure  electron  gas  (such  as  in  gyrotron 
experiments),  we  assume  that  the  electron  gas  is  sufficiently  tenuous 
that  its  static  space-charge  electric  field  can  be  neglected  (a  good 
assumption  for  most  gyrotron  experiments). 

Under  these  assumptions,  the  governing  equations  are  the 
linearized  relativistic  Vlasov  equation, 


3t  fl  + ~ * 3x  fl  " c 3p  fl  = e^l  + c ~x~l^  * 9p  fo  ’ ^ 


9p 


and  the  field  equation 


VxVxE^  = - J Ex  - |t  Jx  , 


(3) 
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where  fQ  and  fj  are,  respectively,  the  equilibrium  and  perturbed 
distributions;  are  the  wave  electric  and  magnetic  fields;  and 

Jl  = - e yVj  vd3p  (4) 

is  the  perturbed  current. 

The  solutions  to  Eqs.  (2)  through  (4)  can  be  obtained  by 
43 

standard  procedures.  Letting  all  quantities  vary  as 
exp(-iu>t  + ik2z)  and  applying  the  method  of  characteristics,  we 
obtain  the  relativistic  dispersion  relation  for  the  right-hand  circularly 
polarized  wave. 


2 u 2 2 

o j - kz  c 


yw  -k2P2/m- 


(5) 


2 2 

where  = 4Trne  /m,  = eBQ/mc,  and 


■ (ltpi 


2/m2c2+pz2/m2ctr2 


f 


(6) 


Equation  (5)  can  be  integrated  by  parts  to  give 


2 


jf  2 .2  2\ 

Pil10  ~ kzc  ) 


272m2c2  (u-k2pz/yfli  - fie/y 


u - k2P2/ym 
oj  - kzpz/ym  - Qg/y 


(7) 
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In  carrying  out  the  integration  by  parts  over  and  pz  the 

integrand  in  Eq.  (5)  has  been  differentiated  with  respect  to  p^  and  pz- 
2 

The  co  term  in  the  numerator  of  the  second  term  on  the  right-hand  side 
of  Eq.  (7)  is  the  net  result  from  differentiations  of  y with  respect  to 
Px  and  pz. 

We  may  now  compare  Eq.  (7)  with  the  standard  nonrel ativistic 

44 

version  of  the  same  dispersion  relation, 


Equation  (7)  is  invariant  in  form  under  Lorentz  transforms  and 

in  the  nonrelativistic  limit  (c— ►oc),  it  reduces  to  Eq.  (7'). 

2 

Comparing  Eqs.  (7)  and  ( 7 1 ) , one  finds  that  the  w term  in 

Eq.  (7),  which  originates  from  the  relativistic  mass  factor  y,  is  absent 

2 2 

from  Eq.  ( 7 ' ) as  expected.  Since  this  term  is  multiplied  by  p^/c  , it 
becomes  small  in  the  nonrelativistic  limit.  However,  it  is  always  an 
important  term  in  the  physical  sense  because  it  gives  rise  to  the  cyclotron 
maser  instability. 

For  the  distribution  function  given  by  Eq.  (1),  Eq.  (7)  reduces 
to 
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where  >0  ^ (l  + pJ/mV)*  and  SiQ  i p^/y^c  . 

For  the  same  distribution  function  in  the  nonrelativistic  limit, 
Eq.  (7' ) reduces  to 


2 i,  2 2 . 

u>  - kzc  = 


k2v  2 
zvlo 


p “ ' Qe  2(u  - flj2 


where  vio  = B1qc 


Again  we  point  out  that  the  difference  between  Eq.  (8)  and 

2 

Eq.  (8 * ) in  the  limit  > 1 i s the  presence  of  the  u term  on  the 

right  hand  side  of  Eq.  (8). 

We  will  now  seek  a unified  physical  interpretation  of  the 
two  instabilities.  Consider  an  electron  moving  with  velocity  vz  in 
the  wave  fields  which  vary  as  exp(-iu>t  + ikzz).  Conventionally,  one 
defines  a Doppler  shifted  wave  frequency  = u>  - kzvz  as  the  effective 
wave  frequency  observed  by  the  moving  electron.  To  facilitate  the 
physical  interpretation  of  the  present  problem,  we  define  instead 
a Doppler  shifted  cyclotron  frequency 


+ 


as  the  effective  cyclotron  frequency  of  a moving  electron  observed 
by  the  propagating  wave. 

Figure  1 shows  the  instantaneous  vector  relationship  of 
the  wave  fields  (JE ^ and  JBj),  the  external  magnetic  field  ( B ) , the 


positions  (points  1 and  2)  and  perpendicular  velocities  (v_^)  of  two 
electrons  located  on  the  same  vertical  line  (x=0)  at  the  initial  time 
t=0.  The  projection  of  the  unperturbed  electron  orbit  on  the  x-y  plane 
is  shown  by  the  dashed  circle.  We  have  assumed  that  the  wave  propagates 
in  the  positive  z-direction. 

The  instantaneous  value  of  at  t=0  is 

%(°)  = k,v2(0)  + ne/y(0)  = ne/y(0).  (10) 

where  we  have  let  vz(0)  = 0 in  accordance  with  Eq.  (1). 

After  an  infinitesimal  time  At,  the  Doppler  shifted  cyclotron 
frequency  fig  will  be 

%(At)  = kzvz(At)  + Se/y(At)  (11) 


or 


Equation  (11)  minus  Eq.  (10)  gives 


Afl0  = kzAvz  + ne  |^orr57  • yXol] 


Afl[j  ~ kzAvz  - (0) 


(12) 


where  AnQ  h nD(At)  - fiQ(0)  , Avz  = v£(At),  diiu  Ay  = y(At)  -y(0). 
From  the  relativistic  equation  of  motion. 


dt 


y v = -eE  - - v x B. 
~ ~ c ~ ~1 


we  obtain 


avz^z-7(^  4x£iAt 


(13) 
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and 


A7  ~ ~ 1 1 1 * £1  At  • 
me 

Substituting  the  relation  (Faraday's  law). 


B,  = — k x E. 

~1  w ~ ~1 

into  Eq.  (13)  to  eliminate  8^,  we  obtain 

Av  = - ( v E. ) At, 

2 y mu  ~L  ~1 
o 

where  we  have  let  y(0)  = yn. 


(14) 


(15) 

(16) 


Substituting  Eqs.  (14)  and  (16)  into  Eq.  (12)  yields 


From  Eq.  (17),  we  draw  the  following  conclusions: 

(i)  Since  v^  for  the  two  electrons  are  oppositely  oriented, 

Eq.  (17)  shows  that  Afig  for  the  two  electrons  will  be  of  opposite  sign. 
Consequently,  if  electron  1 slips  (gains)  in  phase  angle  in  the  phase  space 
of  the  Doppler  shifted  cyclotron  oscillation  (fi^)  electron  2 will  gain  (slip) 
in  phase  angle.  As  a result,  the  two  electrons  tend  to  bunch  toward  each 
other.  The  first  term  on  the  RHS  of  Eq.  (17)  (due  to  Avz)  results  in  axial 
bunching  and  the  second  term  (due  to  Ay)  results  in  azimuthal  bunching. 

(ii)  The  two  bunching  mechanisms  will  be  simultaneously  present 
except  for  the  special  case  k2  = 0 or  = 0. 

(iii)  As  will  be  shown  later,  u>  and  are  of  the  same  sign. 

Hence  the  two  terms  on  the  RHS  of  Eq.  (17)  are  of  opposite  sign,  which 
implies  that  the  two  bunching  mechanisms  always  combine  in  such  a way  as 
to  offset  one  another. 


— ) v,  • Ei 
c2 ) ~-J-  "1 


At 


(17) 


A%  = 


-ek, 


y m 
'o 
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(18) 
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(iv)  The  azimuthal  bunching  mechanism  dominates  if 

con 

— 2T  > 

Yokzc 


I 


and  the  axial  bunching  mechanism  dominates  if 

coO 

bf<*l.  (19) 

Y k c 

T0  2C 

The  same  arguments  leading  to  these  conclusions  also  hold  for 
any  electron  pair  in  the  same  circle  but  on  a vertical  line  with  x f o 
(see  Fig.  1).  Thus,  the  above  conclusions  apply  to  the  system  as  a whole. 

In  comparing  Eqs.  (7)  and  (7'),  we  have  shown  that  even  in  the 
limit  yQ  -*  1,  the  relativistic  and  nonrelati vistic  versions  of  the 
dispersion  relation  of  the  right  hand  circularly  polarized  wave  in  an 
anisotropic  plasma  are  qualitatively  different.  Here  in  the  criterion 
of  Eqs.  (18)  and  (19),  we  again  see  that  o is  not  a decisive  factor 
in  determining  the  relative  importance  of  the  two  bunching  mechanisms. 

This  confirms  our  statement  in  Section  I that  electron  energy  is  not  the 
criterion  to  determine  whether  a relativistic  or  nonrelativistic  model 
should  be  used  to  treat  the  present  problem.  The  proper  criterion  is  qiven 
in  Eqs.  (18)  and  (19). 


Energy  exchange  processes  between  waves  and  hunched  electrons 
for  the  two  instabilities  can  also  be  viewed  in  a unified  way.  Consider 
the  phase  space  of  the  Doppler  shifted  cyclotron  oscillation.  If  the 
electrons  are  uniformly  distributed  in  the  phase  space,  then  clearly  no 
net  energy  can  be  extracted  from  the  electrons  through  their  interactions 
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with  the  wave  electric  field.  If  the  electrons  are  bunched  in  the  phase 
space  of  fig  due  to  either  the  azimuthal  or  axial  bunching  mechanism,  then 
net  beam  energy  loss  is  possible  provided  the  bunched  electrons  remain 
in  the  energy  extracting  phase  of  the  wave  electric  field.  This  imposes 
the  condition  oj  « For  the  distribution  function  of  Eq.  (1),  the 
condition  reduces  to 

w«ne/Y0  (2°) 

Substituting  Eq.  (20)  into  Eqs.  (18)  and  (19),  we  find  that 

fr  > c2  (18') 

Kz 


is  the  regime  where  azimuthal  bunching  dominates,  and 
2 9 

£L  < r2  (io>\ 

kz 

is  the  r lime  where  axial  bunching  dominates. 

alternative  forms,  Eqs.  (18’)  and  ( 19 1 ) , of  the  criterion 
in  Eqs  ) and  (19)  provide  a simple  way  for  mode  identification— the 
azimuth  bunching  mechanism  destabilizes  the  fast  branch  (vacuum  mode) 
of  the  right-hand  circularly  polarized  wave,  while  the  axial  bunching 
mechanism  destabilizes  the  slow  branch  (whistler  mode)  of  the  same  wave. 

One  notes  that,  because  of  the  competitive  nature  of  the  two 
bunching  mechanisms,  only  one  branch  can  be  unstable  for  a given  kz- 
Thus  the  criterion  in  Eqs.  (18’)  and  (19’)  also  serves  to  distinguish 
the  type  of  a given  instability. 
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As  a consistency  check,  we  observe  from  Eq.  (12)  that  for  k2=0, 
only  the  azimuthal  bunching  mechanism  is  present.  Thus  we  expect  an 
instability  from  Eq.  (8)  but  not  from  Eq.  (8‘).  The  absence  of 
instability  in  Eq.  (8’)  for  k2=0  is  obvious.  To  show  that  there  is 
an  instability  in  Eq.  (8)  for  k2=0,  we  neglect  the  first  term  on  the 
right-hand  side  of  Eq.  (8).  The  equation  is  then  readily  solved  to  give 


ne  , . eio“p 

W =S  — + t - r 


This  is  the  cyclotron  maser  instability  in  the  limit  k2*0. 

Substituting  Eq.  (21)  into  the  right  hand  side  of  Eq.  (8), 
we  obtain  the  following  condition  for  neglecting  the  first  term. 


»i.a. 

u\)H  - 


In  the  limit  k2—  °°,  the  axial  bunching  mechanism  dominates 
Tsee  Eq.  (19)].  Then,  both  Eq.  (8)  and  Eq.  (8')  should  yield  the  same 
solution,  as  can  be  easily  demonstrated.  In  this  limit,  u>  is  also  given 
by  Eq.  (21),  except  that  it  is  not  subject  to  the  restriction  Eq.  (22). 


III.  NUMERICAL  EXAMPLES 

In  this  section  we  present  two  representative  numerical  examples 
to  illustrate  the  various  points  made  in  Section  II.  For  reasons 
described  below,  we  may  distinguish  two  characteristic  regimes  according 
to  the  ratio  of  electron  cyclotron  frequency  to  plasma  frequency.  In 
regime  I,  wp  > ^e/yQ  while  in  regime  II,  iop  < ^e/YQ  • 

In  regime  I,  the  cutoff  frequency  (u>p)  of  the  fast  mode 
exceeds  the  electron  cyclotron  frequency.  As  a result,  condition  (20) 
cannot  be  satisfied  for  the  fast  mode  and  only  the  slow  mode  can  be 
unstable  (Weibel-type  instability).  This  is  the  regime  in  which 
nonrelativistic  models  give,  fairly  good  approximations. 

Figures  2a  and  2b  are,  respectively,  solutions  of  Eqs.  (8)  and 

(81 ) (parameters  specified  in  the  caption).  We  have  chosen  a small 

value  for  >0  so  that  the  relativistic  [Eq.  (8)]  and  nonrelativistic 

[Eq.  (8')]  dispersion  relations  are  different  mainly  because  of  the 
2 

presence  of  the  u>  on  the  right  hand  side  of  Eq.  (8).  Comparisons  of 
the  two  figures  show  that  the  nonrelativistic  model,  which  neglects  the 
azimuthal  bunching  mecnanism  in  the  Weibel-type  instability,  only 
slightly  overestimates  the  growth  rate.  In  fact,  the  only  noticeable 
difference  between  the  two  figures  is  the  threshold  k2  for  the 
instability.  The  example  given  in  Fig.  2 is  typical  of  the  entire 
range  of  in  this  regime. 


In  regime  II,  both  types  of  instabilities  exist.  All 
the  conclusions  reached  in  Section  II  are  exhibited  in  this  regime. 
Figures  3a  and  3b  are,  respecti vely,  solutions  of  Eq.  (8)  and  (8') 

l 

(parameters  specified  in  the  caption).  These  two  figures  are  typical  of 
the  entire  T"  demonstrate  the  point  that  relativistic  effect 

may  he  doim.-  • elativistic  energies,  we  have  chosen  a small 

value  for  yQ  (uu  eliding  to  10  keV).  The  dominance  of  relativistic 
effects  results  in  the  destabilization  of  the  fast  mode  (dashed  lines 
in  Fig.  3a).  The  competitive  nature  of  the  two  bunching  mechanisms 
is  apparent  in  Fig.  3a.  We  observe  that  the  growth  rates  drop  sharply 
when  co/kz  « c,  where  the  two  mechanisms  nearly  cancel  each  other.  For  a 
given  k2,  there  is  at  most  one  unstable  branch,  as  has  been  expected 
from  previous  considerations.  Comparing  Fig.  3a  (relativistic  model) 
with  Fig.  3b  (nonrelati vistic  model),  we  find  some  striking  differences 
between  the  two  models.  First,  the  maser  instability  is  absent  from 
the  nonrelativistic  model.  Second,  the  nonrelativistic  model  sub- 
stantially overestimates  the  growth  rate  in  the  region  u>/kz  <,  c , 
where  the  relativistic  bunching  mechanism  is  also  important.  For 
larger  values  of  kz  (io/kz  « c),  the  nonrelativistic  model  begins  to 
yield  more  accurate  results.  Note  that  for  smaller  values  of  kz,  the 
growth  rate  for  the  slow  wave  mode  drops  to  zero  in  the  nonrelativistic 
model.  This  is,  of  course,  not  due  to  cancellation  from  the 
relativistic  bunching  mechanism.  Rather,  it  is  because  condition  (20) 


is  not  satisfied  for  small  values  of  kz< 


IV.  DISCUSSION 


For  the  distribution  function  chosen  [Eq . (1)],  we  have  shown 
that  a nonrelativistic  treatment  of  the  electromagnetic  electron  cyclo- 
tron instabilities  is  essentially  valid  in  the  regime(s)  ^e/Y0wp  < 1 
and/or  u/kz  « c,  but  becomes  very  inaccurate  in  the  regime  ^e/>0ljjp  > 1 
and  ui/kz  < c.  In  the  regime  ®e/y0Up  > 1 and  ui/kz  > c,  nonrelativistic 
models  are  invalid  because  they  completely  miss  the  cyclotron  maser 
instability.  Generalization  of  the  present  results  to  distribution 
functions  with  large  thermal  spread  can  only  be  qualitative  because 
of  the  neglect  of  wave-particle  resonances  in  our  model.  The  principal 
application  of  the  present  theory  appears  to  be  in  the 
area  of  gyrotron  research.  As  noted  earlier,  Eq.  (1)  is  a realistic 
(beam  frame)  representation  of  the  beam  used  in  most  gyrotron  experiments. 
The  physical  relation  between  the  axial  and  azimuthal  hunching 
mechanisms  has  generally  been  overlooked  in  gyrotron  studies.  Clarifica- 
tion of  this  relationship  may  thus  serve  a useful  purpose  in  identifying 
and  predicting  the  type  of  instabilities  in  gyrotron  experiments.  The 
fact  that  only  one  type  of  instability  can  be  present  for  a given  wave- 
length makes  mode  identification  less  ambiguous.  For  example,  the  fast 

4 

wave  generated  in  the  experiment  of  Chou  and  Pantell  has  been  attributed 
to  axial  bunching.  In  the  light  of  the  present  study,  it  appears  that 
the  azimuthal  bunching  mechanism  should  have  been  responsible. 
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Fig.  1 — Instantaneous  (t  = 0)  vector  relationship  of  the  wave  fields  (E1  and  B1), 
the  external  magnetic  field  (B0),  the  position  (points  1 and  2)  and  perpendicular 
velocities  (Vj_)  of  two  electrons.  The  projection  of  the  unperturbed  electron  orbit 
on  the  x-y  plane  is  shown  by  the  dashed  circle,  the  center  of  which  is  taken  to  be 
the  origin  of  the  Cartesian  coordinate  system.  The  positive  z-axis  points  toward 
the  reader. 
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Fig.  2 — (a)  w versus  kz  calculated  from  the  relativistic  dispersion 
relation,  Eq.  (8),  for  70  = 1.02  and  u)e/70c<Jp  = 0.8;  (b)  ui  versus 
kz  calculated  from  the  nonrelativistic  dispersion  relation,  Eq.  (8'), 
for  £2e/wp  * 0.8 


Fig.  3 — (a)  co  versus  kz  calculated  from  Eq.  (8)  for  70  = 1.02 
and  u)e/70a>p  = 10;  (b)  a;  versus  kz  calculated  from  Eq.  (8') 
for  w/cjd  = 10 
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